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Résumé Cet article propose un tour d'horizon des méthodes approchées et exactes, de 
leur performance et de leur complexité théorique, pour différentes versions du problème de 
plus court chemin. L'étude proposée est faite dans l'optique d'améliorer la résolution d'un 
problème plus général de couverture dans le cadre d'un schéma de génération de colonnes, 
dont le plus court chemin apparaît comme le sous-problème. 

Abstract This article provides an overview of the performance and the theoretical com- 
plexity of approximate and exact methods for varions versions of the shortest path problem. 
The proposed study aims to improve the resolution of a more gênerai covering problem 
within a column génération scheme in which the shortest path problem is the sub-problem. 



1 Introduction 

Un problème des plus courants en optimisation combinatoire est celui de la recherche de 
plus courts chemins dans un graphe. Ce problème se présente comme suit : étant donné un 
graphe et une fonction coiit sur les arcs, le problème consiste à trouver le chemin le moins 
coiiteux d'un sommet choisi à un autre. Il se résoud aisément grâce à de nombreux algorithmes 
polynomiaux (Bellman [Bcl58], Dijkstra [Dij59], ...). Néanmoins, l'ajout de contraintes sur le 
chemin (essentiellement, des contraintes de type "sac-à-dos") le rend plus difficile à résoudre. Ce 
rapport de recherche présente donc différents algorithmes exacts ou approchés pour résoudre le 
problème du Plus Court Chemin Contraint (noté PCCC par la suite). Ce travail se place dans 
le cadre plus général de résolution d'un problème de couverture de tâches par des véhicules 
qui doivent respecter certaines contraintes dont le PCCC est, dans une décomposition classique 
de type Dantzig Wolfe, le sous-problème. Aussi, les résultats proposés seront le plus souvent 
comparés dans le cadre du déroulement de tels schémas. Ce document est organisé comme suit : 
la deuxième section présentera le problème PCCC de manière très générale, avec ses différentes 
variantes. Plusieurs méthodes de réduction des instances, qui constituent une étape préliminaire 
à la résolution du problème, sont exposées dans la troisième section. La résolution exacte sera 
abordée dans la quatrième section. Enfin, puisque le problème est NP-difficile, la cinquième 
section portera sur sa résolution approchée. 



2 Description du problème 



2.1 Formalisation 

Soit G = (y, A) un graphe orienté où V est l'ensemble des sommets {\V\ ^ n) et A (- V x V 
l'ensemble des arcs {\A\ ~ m) et soit ^ un ensemble de ressources (|^| = R). A chaque arc 
(i, j) G A sont associés un coût et un vecteur de consommation de ressource (i[j)r=i,...,_R 
positif (Vj' = 1, . . . , iî, t^'j > 0). Le graphe ne doit pas comporter de cycle absorbant"^. 

Un chemin Px^y entre deux sommets x et y est une séquence d'arcs : 
p 

Px-^y = [J{{ut,Vt)} tel que {ui,Vi) G A, Ui+i = Vi, ui = x et Vp = y 
1=1 

La longueur de ce chemin est alors p. Un chemin Px^y allant du sommet x au sommet y a un 
coût dont la formule est C{Px^y) = Cuv 

Soient s et t deux sommets distincts privilégiés du graphe, appelés respectivement source et 
puits. Le problème du Plus Court Chemin Contraint (PCCC) consiste alors à trouver un chemin 
de coût minimal entre la source et le puits satisfaisant certaines contraintes de ressource (ces 
contraintes seront détaillées par la suite). Si l'ensemble 3?. des ressources est vide (iî = 0), on 
est ramené au problème usuel de plus court chemin qui est polynomial ([Bel58], [Dij59]). En 
revanche, la considération d'une seule ressource rend déjà le problème d'optimisation NP-dur, et 
ce même lorsque les coûts et les consommations de ressource sont supposés être entiers positifs 
[GJ79] [Dro94]. Enfin, le problème consistant à décider seulement s'il existe ou non un chemin 
réalisable est NP-complet dès lors que l'on considère deux ressources ou plus. 

2.2 Contraintes de ressource 

Les contraintes de ressource s'expriment à l'aide d'un vecteur consommation de ressource T 
dépendant du chemin considéré dont chaque coordonnée représente une ressource du problème. 

Deux types de contraintes de ressource sont usuellement considérés. Tout d'abord, les contraintes 
de ressource dites finales, où la somme sur tous les arcs du chemin de la source au puits des quanti- 
tés de ressource consommées doit entrer dans une fenêtre définie au puits ; ensuite, les contraintes 
de ressource dites à fenêtres de temps où, à chaque sommet î, sont associées R fenêtres de res- 
source [a\,b\], r = 1, . . . , iî, réduisant l'intervalle des valeurs possibles pour la quantité de res- 
source r pouvant être utilisée avant d'atteindre le sommet i. Pour ce type de contrainte, il existe 
deux méthodes de calcul du vecteur consommation de ressource sur un chemin de la source à un 
sommet : 

- sans attente permise : 

v(F,^, = p,^, u T^{Ps^,) = r{Ps-,ù + tr, 

- avec attente : 

V(P,^, = P,^, U (z, j)), T^{Ps^,) = max{a^^r'-(P,^,) +^1^} 
Un cycle absorbant est un cycle dont la somme des valuations des arcs est négative. 
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Un chemin de la source au puits est dit réalisable si, en chacun de ses sommets, le vecteur 
consommation de ressource sur le sous-chemin de s à ce sommet est dans la fenêtre de ressource 
de ce sommet ; formellement : 

le chemin Ps~tt est réalisable 
ssz V.7 e Ps^t , Vr = 1, . . . , iî, < T^(P,^, ) < 6^ 

De cette inégalité découle deux propriétés générales : 

- dans les deux méthodes (avec ou sans attente) : 

V(ï,j)e A Vr = l,...,iî, <6,^ 

— dans la méthode sans attente permise : 

V(ï,j)e A Vr = l,...,iî, b^ + tl^<a^ 

Les arcs ne vérifiant pas ces inégalités peuvent être directement supprimés car ils n'appartiennent 
à aucune solution réalisable du problème. 

Le premier type de contrainte se ramène aisément au second en attribuant à tous les sommets 
du graphe et pour chaque ressource la fenêtre [0, où 6[ est le majorant de consommation 
définie au puits pour cette ressource. Toutes les instances considérées dans ce document seront 
donc, sans mention expresse du contraire, du second type. 

2.3 Relation avec la génération de colonnes 

Dans certaines modélisations de problèmes linéaires en nombres entiers dites à formulation 
« chemins », les variables du problème représentent des chemins du graphe. Le nombre de chemins 
dans un graphe étant potentiellement exponentiel, il en est de même du nombre de variables. 
L'utilisation d'un algorithme de résolution (simplexe par exemple) ne peut être envisagée pour 
ce modèle puisque ce dernier ne peut être explicité. En revanche, un problème avec un sous- 
ensemble de variables de taille raisonnable, appelé problème maître restreint (PMR), peut être 
résolu ; la résolution de ce dernier permet de calculer les coiits réduits associés à chaque arc 
du graphe, que l'on peut interpréter comme le coût d'opportunité à emprunter un arc donné 
(par le biais d'un chemin utilisant cet arc). La variable pouvant être ajoutée au PMR est alors 
trouvée par résolution du sous-problème qui consiste en la détermination d'un plus court chemin 
contraint pour le critère de coiit réduit. Par exemple, dans le cas d'une minimisation, si le coût 
réduit total d'un chemin est négatif, il est potentiellement améliorant pour le problème maître et 
constitue ainsi un bon candidat à être intégré au PMR. Dans cette approche appelée génération 
de colonnes, l'optimalité est atteinte lorsqu'il n'existe plus de chemin améliorant (chemin de coût 
réduit négatif si l'on considère la résolution continue du problème maître). Notons que, pour la 
convergence, il n'est pas nécessaire de trouver un chemin optimal : un chemin améliorant suffit. 
De plus, à chaque itération, plusieurs chemins améliorants peuvent venir enrichir le PMR. 
Le problème de plus court chemin contraint étant NP-difïicile, le temps de résolution des sous- 
problèmes n'est pas maîtrisé ; d'où l'intérêt d'utiliser des algorithmes d'approximation, du moins 
en début de schéma, nous permettant de trouver rapidement des chemins améliorants (la rsolution 
exacte demeurant inévitable en fin de schéma pour prouver l'optimalité). 

3 Réduction du problème avant sa résolution 

Quelques travaux précédant la résolution peuvent être effectués dans le but de réduire le 
graphe (suppression de sommets ou d'arcs, réduction de l'amplitude des fenêtres de temps). En 
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outre, ces traitements permettent parfois de détecter les instances non réalisables {i.e., instances 
sur lesquelles tout chemin viole les contraintes de ressource) et d'exhiber un majorant et un 
minorant. 



3.1 Réduction des fenêtres de temps 



Les deux types de contraintes de ressource forment en fait une seule classe de problèmes oii 
chaque sommet possède des fenêtres de temps. En chaque sommet, les bornes des fenêtres de 
temps pour chaque ressource doivent vérifier quelques relations dépendant de ses prédécesseurs 
et de ses successeurs. L'ensemble des prédécesseurs de i se notera pred{i), l'ensemble des ses 
successeurs succ(i). Ainsi, pour une ressource donnée r, le minorant d'une fenêtre du sommet j 
ne doit pas être plus petit que le minorant d'un prédécesseur i du sommet auquel le temps de 
trajet entre les sommets i et j est ajouté. Cela donne l'équation suivante : 



= max{a^, min {a[ + 

iepredÇj) 



Une équation similaire est valable pour les majorants : 

6^ = min{6^, max + 



Maintenant, en considérant les successeurs, le minorant d'un sommet j ne doit pas être plus petit 
que le minorant d'un successeur i du sommet auquel le temps de trajet de j à i est soustrait. 
Cela donne l'équation suivante : 



= max{a^, min {a\ — t^J} 

iesucc{j) 



Une nouvelle équation similaire est valable pour les majorants : 



b'j = min{ 6;' , max { fe[ - t } } 



Si, pour un sommet i et une ressource r, la condition > b\ est remplie, alors le sommet i doit 
être supprimé du graphe (car il est inaccessible), ainsi que tous ses arcs incidents. De même, dans 
le cas où l'attente n'est pas permise, s'il existe un arc (i, j) tel que b\ + < aj, alors l'arc (i, j) 
peut être supprimé. 



L'algorithme 1 vérifie donc les équations précédentes pour chaque sommet. Dans le cas d'un 
graphe acyclique, il suffit de traiter les sommets dans l'ordre topologique pour les équations utili- 
sant les prédécesseurs, puis dans l'ordre inverse pour les équations concernant les successeurs. La 
procédure recommence si et seulement si un sommet ou un arc peut être supprimé. La complexité 
de cette procédure est donc en 0['n? R). Dans le cas cyclique, rien ne garantit qu'aucune fenêtre 
ne pourra encore changer après l'application de la procédure, et ce même si celle-ci n'induit pas de 
suppression de sommets ou d'arcs. La procédure précédente est donc répétée tant qu'une borne 
change ou qu'un sommet ou un arc est éliminé. Pour cet algorithme, la complexité en temps 
de la boucle répéter est en 0{v? R). Une première approximation naïve du nombre de fois oit 
cette boucle est effectuée dans le cas de fenêtres entières est 0{n max{6i — a;}) (réduction d'au 

moins une unité de la fenêtre pour un sommet à chaque itération) . La complexité totale est donc 
0{n^ R max{6i - aj). 
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Algorithme 1 : REDUCTION : Réduction des fenêtres de temps (graplie acyclique) 
Entrées : 

G = {V, A) et Vi € V^, Vr = 1, . . . , R, [al; bl] la fenêtre de la ressource r pour le sommet i ; les sommets 
sont numérotés dans l'ordre topologique 

Sorties : 

G' — {V',A') un sous-graphe partiel de G avec les fenêtres de ressource réduites 

V ^V; A' ^ A; 
répéter 

recommencer ^ false; 
pour i allant de 1 à n faire 

pour r allant de 1 à R faire 
minp ^ oo; maxp <— 0; 
pour tous les j £ pred{i) faire 
minp ^ min {minp, + t^ji}', 
maxp ^ max{maxp, + tj^}; 

si a[ < minp alors al ^ minp; 
si bl > maxp alors bl <— maxp; 

pour i allant de n à 1 faire 

pour r allant de 1 à R faire 
mins ^ oo; maxs ^ 0; 
pour tous les j £ succ{i) faire 
mins ^ raxa{mins , aj — tj^}; 
maxs <— miLK{maxs, bj — tj^}; 
si al < mins alors al <— mins; 
si bl > maxs alors bl ^ maxs; 

si attente non permise et 3(i,i) G A' , 3r = 1, . . . , R,bl + tlj < Oj alors 

|_ recommencer <— true; 

si 3i £ V' , 3r = 1, . . . , R,bl < al alors 
V'<-V'\{^}; 

A' ^ A' \ {{u,v) £ A' ,u ^ i ou V ^ i}; 
recommencer <— true; 

jusqu'à recommencer — false; 
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3.2 Prétraitement pour contraintes de ressource finales 

Le but de cette procédure de prétraitement, présentée dans [AAN83] et améliorée dans[DB03] 
est, comme précédemment, de réduire le graphe mais en plus de fournir un minorant et un ma- 
jorant pour la valeur du problème. Elle exploite le fait suivant : trouver le plus court chemin non 
contraint entre deux sommets a la même complexité (et prend aussi le même temps de calcul) 
que trouver le plus court chemin entre un sommet et tous les autres (ou tous les autres et un 
sommet). Les plus courts chemins suivant chaque métrique (coût ou ressource) sont donc cal- 
culés entre la source et tous les sommets, ainsi qu'entre tous les sommets et le puits. Ensuite, 
en recombinant ces plus courts chemins, certains peu coiiteux mais non réalisables peuvent être 
détectés, ce qui permet d'améliorer le minorant ; inversement, des solutions réalisables peuvent 
être exhibées, ce qui fournit un majorant au problème. Enfin, ce traitement permet d'élaguer le 
graphe en déterminant des arcs et des sommets qui n'appartiennent à aucun chemin réalisable 
ou qui n'appartiennent à aucun chemin optimal. 

Le déroulement de l'algorithme est le suivant : 

- Tout d'abord, les chemins les moins coûteux de la source à tous les sommets sont calculés. 
Cela permet éventuellement de détecter des instances non réalisables (absence de chemin 
de la source au puits) ou d'exhiber un chemin optimal (un plus court chemin de la source 
au puits vérifie les contraintes de ressource) ou encore, d'exhiber un chemin non réalisable 
de coût minimum permettant de mettre à jour le minorant pour le problème. 

- Ensuite, pour chaque ressource, les chemins les moins consommateurs de cette ressource 
allant de la source à chaque sommet sont calculés. Ils permettent une nouvelle fois de détec- 
ter la non faisabilité de l'instance (la valeur du plus court chemin sur une ressource dépasse 
le majorant du puits) ou de mettre à jour le majorant du problème par la considération 
d'un chemin réahsable pour toutes les ressources. 

- Enfin, les chemins les moins coûteux et les moins consommateurs en chaque ressource allant 
de tout sommet au puits sont calculés. 

D'une part, l'ensemble des chemins ainsi générés permet d'améliorer le majorant en recombinant 
les plus courts chemins pour trouver des chemins réalisables de bon coût : le principe consiste 
tout simplement à considérer, pour chaque arc (î,j) du graphe, les chemins résultant d'une 
concaténation d'un plus court chemin de s à i, de l'arc (i, j) et d'un plus court chemin de j à t. 

D'autre part, cela permet de tester l'accessibilité ou la pertinence (en incluant l'information 
apportée par le majorant) de chaque sommet et de chaque arc pour les supprimer si possible : 
pour un sommet i {resp., pour un arc il suffit que, pour une métrique donnée, la valeur 

du plus court chemin de s à i additionné à la valeur du plus court chemin de i à t {resp., à la 
valeur de l'arc {i, j) plus la valeur du plus court chemin de j à t) soit plus grande que le majorant 
en coût ou que le majorant présent sur le puits pour cette métrique pour que ce sommet (resp., 
cet arc) soit supprimé. Cette procédure peut donc renvoyer un constat de non réalisabilitc, un 
chemin optimal ou un minorant et un majorant, ce dernier pouvant être associé à un chemin. 
Une description de ce traitement est proposée par l'algorithme 2, où l'on suppose disposer de 
deux procédures pcc{s — >, /) et pcc{-^ t, /) qui permettent de déterminer respectivement les 
plus courts chemins de s à tous les sommets et les plus courts chemins de tous les sommets à t, 
relativement à la métrique /. 

Dans le cadre d'un schéma de génération de colonnes, ce traitement permet de détecter 
une instance non réalisable mais aussi d'éliminer des sommets et des arcs n'appartenant pas à 
une solution optimale. De plus, dès que la valeur du majorant devient négative, l'élimination de 
sommets et d'arcs non optimaux n'est plus nécessairement opportune si l'objectif est de conserver 
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plusieurs solutions intéressantes (de coût négatif). Il suffit alors de rendre stricts les tests sur le 
majorant dans les blocs d'élimination des sommets et des arcs. 

4 Résolution exacte du PCCC 

Nous rappelons que le problème PCCC est NP-dur, et ce même pour une seule ressource. Le 
problème de décision associé à l'existence d'un chemin de coût inférieur à une borne est lui-même 
NP-complet pour le cas de deux ressources ou plus. La résolution exacte de ce problème peut 
être menée par la programmation dynamique. 

4.1 Cas général 

La programmation dynamique pour ce problème permet d'élaborer un algorithme pseudo- 
polynomial^ de résolution exacte. La programmation dynamique se fonde sur le principe d'op- 
timalité de Bellman : toute sous-séquence d'une séquence optimale est optimale. Cela est di- 
rectement appliquable au plus court chemin : si P = {(s,vi), . . . , {vq,t)} est un plus court 
chemin de s à alors P' = {(s, vi), . . . , {vi-i, Vi)} est nécessairement un plus court chemin de 
s h Vi pour i < q. Avec l'introduction des contraintes de ressource, on ne peut plus adapter 
directement ce principe pour propager les meilleurs chemins, puisqu'il faut nuancer la notion 
d'optimalité des sous-séqucnccs par la considération d'un niveau de consommation de ressource : 
si P ~ {(s, wi), . . . , {vq, t)} est un plus court chemin de s à t consommant {B^ , . . . , B^) unités 
des R ressources, alors P' = {(s,wi), . . . , (vi-i,Vi)} est un plus court chemin de s à Vi, parmi 
les chemins ne consommant pas plus de {B^ — T^{Pi^t), ■ ■ ■ ,B^ — T^{Pi-,t)) unités des res- 
sources. Pour gérer les niveaux de consommation de coût et de ressource des sous-séquences, la 
programmation dynamique utilise la notion d'étiquettes. 

Définition 1 Étiquette 

Une étiquette est un vecteur représentant un chemin et dont les coordonnées sont le coût et les 
consommations des différentes ressources. 

E={E^,E'^,..., E^) où E° est le coût du chemin et E'' pourr G {1, . • . , R} est la consommation 
de la ressource r. 

Remarque 1 

A tout chemin réalisable entre s et tout sommet correspond une étiquette. 

Les différents algorithmes de programmation dynamique utilisent ces étiquettes, mais il n'est pas 
forcément nécessaire de toutes les garder : conserver toutes les étiquettes reviendrait à énumérer 
tous les chemins du graphe. Pour éliminer les étiquettes inutiles, une relation de dominance est 
définie. 

Définition 2 Relation de dominance : 

Une étiquette E domine une étiquette E' ssi Vr € {0, . . . , R}, E^ < E'^ 
{E >- E') et 3r e {0, . . . , R}, E"- < E"^ 

Cette relation de dominance induit un ordre partiel appelé ordre de Pareto. Toutes les étiquettes 
n'étant pas forcément deux à deux comparables, cela permet de définir un ensemble de majorants 
pour cet ordre. 

* La complexité dépend polynomialement de la taille des instances ainsi que des données numériques. 
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Algorithme 2 : PRETRAITEMENT : Procédure de prétraitement 
Entrées : 

- G = {V, A) graplie, c : A ^ R fonction de coût sur les arcs, b G vecteur des ressources disponibles 
au puits 

Sorties : G' sous-grapiie partiel de G ; U et L majorant et minorant du chemin le moins coûteux et 
vérifiant les contraintes 

V ^V; A' ^ A; 0; 

f/ <— f/o = Cmax X (1^1 - I) + 1 avec Cmax = max {Cij}; 

(iJ)6A 

répéter 

chg false; 

Il Calcul des plus courts chemins 

{Ps^i.i G V"} ^ pcc(s ^,c); 

Vr = 1, . . . {PU,,i € V'} ^ pccis 

Vr = 1, . . . , iî, {p;L,„i e V'} ^ pcc{^ t, n- 

Il Considération de Ps—>t 

si $Ps^t alors L'instance n'est pas réalisable ; EXIT; 
sinon si Vr = 1, . . . , iî, T''[Ps^t) < bl alors 
|_ Ps~tt est une solution optimale ; EXIT; 

sinon L ^ C{Ps~,t); 
Il Considération de 

si 3r, r''(PJ_() > 6[ alors L'instance n'est pas réalisable ; EXIT; 

si 3r, = 1, . . . , R, T%P:^,) < 6? et C(P;_0 < U alors U ^ C{P^^t); 

Il Recombinaison des chemins 
pour tous les G A' faire 
si 

R R 

3PG{A^.}uy{PiL.a, 3P' e{p;^t}u\j{p;r,,}, w = i,...,r, r\p)+t:,+T'-{p')<bi 

|_ et G{P) + + C(P') < U alors U ^ C{P) + + C(P'); 

// Élimination de sommets 
pour tous les i £ V' \ {s, t} faire 

si 3r = 1, . . . , P, T'-(PJ_,) + T^iPr^t) > bl alors 
|_ Supprimer le sommet i et tous ses arcs incidents; chg ^ true; 

si C{Ps-,i) + C{Pl^t) > U alors 
|_ Supprimer le sommet i et tous ses arcs incidents; chg ^ true; 

Il Élimination d'arcs 
pour tous les £ A' faire 

si 3r = 1, . . . , P, T'-(P;_^,) + tl^ + r'-(P;^,) > 61 alors 
|_ Supprimer l'arc (i, j); chg ^ true; 

sinon si C{Ps-,i) + c,j + C{P'j_t) > U alors 
|_ Supprimer l'arc (i, j); chg <— true; 

jusqu'à chg — false; 

Retourner les bornes L et U ainsi que le chemin correspondant à la borne U; 
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Définition 3 Élément Pareto- optimal ou non dominé : 
E est Pareto- optimal ssi $E' , E' >- E 

Définition 4 Ensemble Pareto- optimal : 

Un ensemble Pareto- optimal est un ensemble d'éléments non dominés. 

Pareto = {E \ $E' , E' >- E} 

Cette relation de dominance permet de ne générer que les chemins Pareto-optimaux. Pour fc € N, 
la conservation de toutes les étiquettes non dominées par k autres permet de s'assurer de trouver 
les k meilleures solutions et permet ainsi, dans notre schéma de génération de colonnes, d'insérer 
plusieurs colonnes lors d'une même itération. 

Deux types d'algorithmes de programmation dynamique existent : algorithmes à correction 
d'étiquettes ([DPS83]) et algorithmes à fixation d'étiquettes ([DS88]). 

Algorithme à correction d'étiquettes [DPS83] 

Cet algorithme se nomme ainsi car, à chaque itération, il va essayer d'améliorer les étiquettes 
déjà existantes. En effet, une liste de sommets sur lesquels il existe des étiquettes non encore 
traitées est maintenue. A chaque étape, un sommet de cette liste est choisi et les étiquettes de ce 
sommet sont propagées aux successeurs de ce sommet. Pour chacune des étiquettes ainsi créées, 
un test de dominance est effectué pour éliminer les étiquettes dominées ou déjà obtenues. Si une 
nouvelle étiquette apparaît, le sommet sur lequel cette étiquette est présente est ajouté à la liste 
des sommets à traiter. 



Algorithme 3 : CORRECTION : Algoritlime à correction d'étiquettes 
Entrées : 

- G = {V, A) : graphe avec R contraintes de ressource 

- Pareto{S) renvoie l'ensemble des étiquettes non dominés de S 

Sorties : Chemin contraint de coût optimal 

// ETIQ{i) est l'ensemble des étiquettes du sommet i, les étiquettes étant des vecteurs 

de taille (1 + R) 
pour tous les i £ V faire 

L ETIQ{i) ^ 0; 

LIST ^ {s}; 

ETIQ{s) ^ {0} 

tant que LIST / faire 

Choisir i G LIST; LIST ^ LIST \ {i}; 
pour tous les j G suœ{i) faire 

pour tous les E G ETIQ{i) faire 

si Vr G {1, . . . , -R}, + tlj < alors 

E' ^{E° +c,,,E^ +t1^,^re{l,...,R}); 
ETIQij) <- Pareto{ETIQ{j) U {£'}); 
\_ si E' e ETIQ(j) alors LIST ^ LIST U {j}; 

Retourner le chemin ayant le plus petit coût en t; 
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Algorithme à fixation d'étiquettes [DPS83] 

Cet algorithme se nomme ainsi car à chaque itération, il va fixer une étiquette qui ne pourra plus 
être modifiée. En effet, une liste d'étiquettes non encore traitées est maintenue. A chaque étape, 
une étiquette (une des non dominées de la liste) est choisie pour être propagée. Les nouvelles 
étiquettes créées, si elles ne sont pas dominées sur leur sommet, sont rajoutées à la liste. 



Algorithme 4 : FIXATION : Algorithme à fixation d'étiquettes 
Entrées : 

- G = {V, A) : graphe avec R contraintes de ressource 

- Pareto{S) renvoie l'ensemble des étiquettes non dominés de S 

- Sommet{E) renvoie le sommet d'étiquette E 

- min_ordreJex{S') renvoie un élément de l'ensemble <f minimum pour l'ordre lexicographique 

Sorties : Chemin contraint de coût optimal 

// Les étiquettes sont des vecteurs de taille (1 + iî). 
// ETIQ(i) est l'ensemble des étiquettes du sommet i. 

Il DEF_ETIQ{i) est l'ensemble des étiquettes définitives du sommet i. 
pour tous les i e V faire ETIQ{i) ^ 0; DEF_ETIQ{i) ^ 0; 

ETIQ(s)^{{0,...,0)} 

tant que (J {ETIQ{i) \ DEF_ETIQ{i)) / faire 

iev 

E ^ min_ordreJex{ [J {ETIQ{i) \ DEF_ETIQ{i))); 

i «— sommet{E); 

pour tous les j £ succ(i) faire 

si V£ G {1, . . . , iî}, -E* + 4j < h] alors 

E' ^ {E^ + Cij ,E^ +tjj, . . . , E^ +tfj); 
|_ ETIQij) ^ Pareto{ETIQ{j) U {£'}); 

_ DEF_ETIQ{i) ^ DEF_ETIQ{i) U {E}; 
Retourner le chemin ayant le plus petit coût en t; 



4.2 Cas d'un graphe acyclique 



Dans le cas d'un graphe acyclique, un ordre topologique sur les sommets peut être calculé. 
L'ensemble des sommets peut alors être numéroté de 1 à n de sorte que le sommet 1 soit la 
source, le sommet n le puits et que les sommets de tout arc vérifie la relation i < j. Cet 
algorithme de programmation dynamique pour les graphes acycliques va appliquer l'algorithme 
à correction d'étiquettes vu précédemment en parcourant les sommets dans l'ordre topologique. 
Ainsi, lors du traitement d'un sommet, tous les sommets qui le précèdent dans l'ordre topologique 
ont leurs étiquettes définitives. Le calcul des étiquettes de ce sommet ne se fait donc qu'une seule 
fois. 



Etude de complexité 

Pour calculer la complexité en espace, il suffit de calculer le nombre maximum d'étiquettes que 
peut générer cet algorithme. En supposant cpie chaque sommet est relié à tous les sommets le 
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Algorithme 5 : ACYCLIQUE : Algorithme pour les graphes acycliques 
Entrées : 

- G — {V, A) : graphe acyclique avec R contraintes de ressource, V — {1 . . .n} 

- Pareto{S) renvoie l'ensemble des étiquettes non dominés de S 

Sorties : Chemin contraint de coiit optimal 

// ETIQ{i) est l'ensemble des étiquettes du sommet i, les étiquettes sont des vecteurs 

de taille 1 + R 
ETIQ(1)^{{0,...,0)}; 
pour tous les i allant de 2 à n faire 
ETIQ{i) ^ 0; 

pour tous les j e pred{i) faire 

pour tous les E £ ETIQ{j) faire 

si e {1, . . . , R}, E^ + t% < b^j alors 

E' ^ {E" +c,,,E' +t}^,...,E'' + tfjy, 
ETIQ{i) ^ ETIQ{i)U {E'}); 

ETIQ{i) Pareto{ETIQ{i))- 
Retourner le chemin ayant le plus petit coiit au sommet n; 



précédant dans l'ordre topologiquc, le nombre maximum d'étiquettes pour le sommet i (Ni) est 
égal au nombre de toutes les étiquettes des sommets de plus petit indice ; la source n'ayant 
qu'une seule étiquette (nulle) : 



■i-i 

Ni = let N,, ^^Nj, i>l 
soit, Ni^let N, = 2*-^ i > 1 



Pour trouver le nombre d'étiquettes total, il suffit de faire la somme sur tous les sommets 



i=l i=2 



Sachant que chaque étiquette est un vecteur de taille 1 + iî, la complexité en espace est donc en 
0{R e"). 

La complexité en temps de l'algorithme précédent se calcule en prenant le cas du sommet 2 
à part car il n'y a qu'une seule étiquette créée et donc pas de dominance (i itère les sommets, j 
les prédécesseurs, k les étiquettes, i et m sont deux itérateurs pour faire les comparaisons deux 
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à deux lors de la dominance) : 



1=2 \7=1 




R + N2R + C^.iî + R^\Y,Nj+G\ 



i=3 y=i 

2^~^ j < 2, on obtient : 



n / i—\ 



en sachant que A^i = 1, A^2 = 1 et N.j 



i=3 ^ ' 




+ 



Ni-l Ni 



E E 



R 



e") 



î=3 \i=i fe=i 



L'algorithme de programmation dynamique proposé pour les graphes acycliques n'est pas 
polynomial. Les algorithmes pour les graphes généraux ne le sont a fortiori pas non plus. Au vu 
de ces résultats de complexité, l'élagage du graphe a priori en supprimant des sommets ou des arcs 
est primordial pour réduire le nombre de chemins possibles ; de même, diminuer l'amplitude des 
fenêtres de temps et trouver un bon majorant permet aussi de diminuer le nombre d'étiquettes 
calculées. 

4.3 k Plus Courts Chemins 

Des algorithmes polynomiaux exacts ont également été développés pour déterminer, non plus 
un plus court chemin, mais les k plus courts chemins de s à i (ou de s à tout autre sommet). 
Le meilleur algorithme connu à ce jour est celui de Eppstein, [Epp98], de complexité 0{m + 
n log(n) + k). Il consiste à fabriquer un tas des chemins du graphe stockés sous forme implicite 
et trié selon un certain critère. L'explicitation d'un chemin du tas se fait ensuite en 0{p) on p 
est la longueur du chemin à énumérer. 

Dans le cas du PCCC, cet algorithme peut calculer les consommations de ressource s'addi- 
tionnant le long de tous les chemins lors de la construction du tas avec une complexité d'ordre 
0{m + Rn log(n)). La récupération de la valeur d'une consommations se fait en temps constant. 
Ainsi, il est possible de tester la réalisabilité d'un chemin en temps 0{R). 

Pour résoudre le problème par cet algorithme, il faudrait générer le tas de tous les chemins 
du graphe selon le critère de coût initial, puis dépiler les chemins (dont le nombre total est 
exponentiel) jusqu'à trouver un chemin réalisable, le premier chemin trouvé étant le chemin 
optimal. Une telle procédure, de complexité au pire des cas d'ordre exponentiel en temps, pourrait 
néanmoins consituter une alternative à la résolution exacte par programmation dynamique. 

5 Résolution approchée du PCCC 

Le problème du PCCC étant NP-difRcile, la recherche s'est dirigée vers la résolution ap- 
prochée du problème. Quelques heuristiques existent : aggrégation de contraintes de ressource 
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[NS05], recherche de solutions e-réahsable [ABS02]. Cette section présentera un schéma d'ap- 
proximation pour le cas d'une ressource finale puis une généralisation de ce schéma dans le cas 
de plusieurs ressources. 

5.1 Schéma d'approximation totalement polynomial 

Un FPTAS [Fully Polynomial Time Approximation Scheme ou schéma d'approximation tota- 
lement polynomial en temps) a tout d'abord été présenté par Hassin [Has92] et ensuite amélioré 
par différents auteurs [Phi93,LR01;ESZ02]. Ce FPTAS se base sur la programmation dynamique 
avec la technique d'échelonnage et d'arrondis pour trouver une approximation de la solution. Ce 
FPTAS est utilisable dans le cas de graphes acycliques, avec une seule contrainte de ressource 
finale. Un ordre topologique sur les sommets est calculé ; l'ensemble des sommets est donc doré- 
navant numéroté de 1 à n sachant que le sommet 1 est la source, le sommet n le puits et tout 
arc (î, j) vérifie la relation i < j. 

Schéma initial La procédure de programmation dynamique sous-jacente est basée sur le coût 
et non sur la ressource, comme c'est le cas habituellement : on cherche à déterminer pour c = 
1, 2, . . . , gj(c) la plus petite consommation de ressource des chemins allant de 1 à j en au plus c 
unités de coiit ; la recherche porte alors sur gn{c), c étant optimal des lors que 5n(c) < (i-C, 
la valeur optimale est le plus petit c pour lequel la consommation de ressource est inférieure à la 
borne du puits). L'utilisation de cette procédure induit donc que les coiits sur les arcs sont des 
entiers strictement positifs. 



Algorithme 6 : EXACT : Programmation dynamique basée sur le coût 
Entrées : Un graphe G acyclique 

Sorties : Un chemin contraint de coût optimal (OPT) 
0; 

51 (c) ^0; 

pour tous les j £ {2, . . . , n} faire gj{c) ^ oo; 
tant que gn(c) > b„ faire 
c+1; 

pour tous les j G {2, . . . , n} faire 

5j(c) ^ gj{c- 1); 
pour tous les G A faire 

|_ si c,j < c alors gj{c) <- mïn {gj{c), gi{c - dj) +tij}; 

Retourner le chemin de coût c correspondant à la consommation de ressource (?n(c) 



La complexité de l'algorithme EXACT est 0{m OPT) avec n ^ m. Il s'agit donc d'un 
algorithme pscudopolynomial, i.c., dont la complexité est polynomiale en la taille de l'instance 
prise au sens du nombre d'éléments de la structure à coder, mais exponentielle en le logarithme 
de ses données numériques. 

Une manière usuelle de se ramener à un ordre polynomial de complexité consiste à diminuer 
considérablement l'ordre de grandeur des données numériques ; c'est la technique d'échelonnage 
et d'arrondis. Bien sûr, la polynomialité de l'algorithme appliqué à l'instance transformée se 
gagne au prix de l'optimalité : les solutions obtenues ne sont plus optimales, mais seulement 
approchées, pour l'instance initiale. 
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Pour le problème qui nous concerne, on utilise l'algorithme SCALING qui, étant donnés une 



borne £? G N et un rationnel ô g]0, n], remplace chaque coiit Cy par le coiit 



(n- 1) 
ô B 



Algorithme 7 : SCALING ; Procédure d'échelonnage et d'arrondis 
Entrées : I = (G = {V, A), c : A ^ N) graphe arc-valué, B eN, â g]0, n] 
Sorties : I{B,S) graphe arc-valué 

0; 

pour tous les £ A faire 
si Cij < B alors 

Cîj(n- 1) 

L SB y 

|_ À^Àu{i,j); 



La complexité de la procédure SCALING est en O (^m log (^'g^^ (pour chaque arc, une re- 

cherche dichotomique dans l'espace [0, — ]). De plus, la procédure EXACT appliquée à l'instance 



I { Ij. ô) nécessiterait un temps 0{m OPT), où OPT désigne la valeur optimale sur I{B, ô). Or, les 
valeurs de toute solution Pg^t sur les instances / et I{B, ô) sont liées par les relations suivantes : 

CiPs^t) < -^C(Ps^t) + Bô C{Ps-.t) < ^C'iPs^t) 

n — l Bà 

En particulier, si les chemins Ps—,t et P*^^ sont respectivement optimaux pour les instances 
I{B, S) et /, alors : 

~ n — 1 

OPT < C(P*f) < -OPT 

Bô ^—^^ 

C(Ps^t) < OPT + BÔ < OPT + Bô 

n — 1 

On déduit de ces relations : d'une part, si OPT est à rapport polynomial p de B (i.e., 
OPT < pB où p est borné par un polynôme en n), alors le déroulement de l'algorithme EXACT 
sur I{B,ô) devient polynomial; d'autre part, si B est un minorant de OPT (i.e., B < OPT), la 
solution renvoyée est (1 + (5)-approchée pour le problème initial. Il revient donc à déterminer en 
temps polynomial une borne B vérifiant B < OPT < pB. 

La détermination de B s'effectue par recherche dichotomique dans l'intervalle [LB, UB], où 
LB et UB désignent respectivement un minorant et un majorant de OPT (par exemple, considé- 
rer LB = 1 et UB ~ (n— 1) X max {cu}). Cette recherche est elle-même fondée sur la procédure 

TEST qui, pour une instance / et une borne B données, renvoie : 



TEST{B,Ô) = 



OUI si OPT > B 
NON si OPT < (1 + ô)B 



Une mise en œuvre de cette procédure consiste à appliquer la programmation dynamique EXACT 

7j — 1 

sur l'instance I{B, ô), en limitant le coût maximum à tester à — - — . 
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Algorithme 8 : TEST : Procédure de test approchée 



Entrées : / instance, B borne sur le coiit objectif, ô €]0, n] erreur 
Sorties : OUI si OPT > B ou cliemin P»-.* de valeur C(Pa^t) < (1 + S)B 
Î{B,Ô) ^SCALING(/,B,(Î); 




pour tous les c de 1 à B faire (;s(c) ^ 0; 
pour tous les j £ {2, . . . , n} faire gj{0) <— oo; 
pour c de 1 à B faire 

pour tous les j G {2, . . . , n} faire 

^ ffj(c- 1); 
pour tous les G A faire 

|_ si Cjj < c alors ^^(c) ^ min {gj (c), pi(c - c^j) 

si g,i(c) < bn alors Retourner NON et le chemin trouvé; 
Retourner OUI 



La complexité de la procédure TEST est dominée par l'étape de programmation dynamique, 

Tïl Tl 

qui s'effectue en temps 0{ — —). Lorsque la procédure TEST répond NON, le chemin renvoyé, 



de valeur au plus ^ — - — ^ dans l'instance Î{B, S), est de valeur au plus B{1 + ô) dans le graphe 

initial. Si, en revanche, elle répond OUI, c'est que tout chemin réalisable est de valeur au moins 
n — l ~ 

— : h 1 sur I(B, ô), et donc, de valeur au moins B sur /. Autrement dit, la décision exacte sur 

à 

I se transforme en décision (5-approchée sur /. 

~ i OUI ^ OPT > è ( OPT > B 

\ NON ^ OPT <É^ l OPT < (1 + S)B 

La recherche dichotomique DICHO consiste alors, partant d'un intervalle initial [LB, UB], à 
appeler itérativement la procédure TEST, jusqu'à obtenir un encadrement suffisamment fin de la 
valeur de OPT {UB < pLB). Si l'on pose f{LB, UB) = (5 G ]0, n] et g{LB, UB, ô) = ^JLB x UB, 
la complexité de la procédure DICHO est en 

oii O l log ( — tA ^ I I estime le nombre de tests nécessaires (recherche dichotomique entre 

V V log(p) // 

LB et UB dans l'espace logarithmique) et C'(loglog(J7i?/Li3)) le temps de calcul d'une valeur 
approchée (mais suffisante) de g{LB, UB, S). Cette valeur approchée est trouvée en cherchant le 

premier indice i tel que 2^ > - ; on prend alors g{LB, UB,5) = LB 2^ 

LB 

Le schéma d'approximation totalement polynomial peut maintenant être énoncé. Pour une 
instance / et une erreur e e]0, 1[ données, le schéma se décompose en trois étapes : 

- Tout d'abord, il calcule un minorant LB et un majorant [/S de la valeur du problème. 

- Ensuite, il raffine à l'aide de DICHO l'encadrement de l'optimum jusq'à optenir LB < 
OPT <UB < pLB, p e K 

- Enfin, il détermine le chemin optimal dans l'instance I{LB, e) grâce à EXACT. 
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Algorithme 9 : DICHO : Recherche de la borne B 



Entrées : LB et UB bornes inférieure et supérieure de OPT, p paramètre 
Sorties : Une borne LB telle que LB < OPT < pLB 

tant que UB > pLB faire 

5 f{LB,UB); 
B ^ g{LB,UB,S); 

si TEST{B,S) répond OUI alors LB ^ B; 
_ sinon UB ^ B{1 + 5); 



Algorithme 10 : SCHEMA : Schéma d'approximation totalement polynomial 
Entrées : e l'erreur d'approximation, p paramètre 
Sorties : Un plus court chemin contraint (1 + e)-approché 

(1) Déterminer LB et UB bornes inférieure et supérieure de OPT; 

(2) LB ^mC}10{LB,UB,p); 

(3) Ps^t ^EXACT(/(LB,£)); 



La complexité des étapes (2) et (3) de cet algorithme est donc de 

K'-(^^)(^--'-(fi)))-("'H)) 

Hassin [Has92] utilise les paramètres constants lî = e et p = 2 pour la procédure de réduction 
de l'intervalle [LB,UB] ; il initialisc par ailleurs les bornes LB et UB respectivement à 1 et 
{n — l)Cmax- La complexité finale du schéma qu'il propose est ainsi d'ordre : 

-(>--(fl)(^ 

qui est bien polynomiale en (n, m, log(Cmaa:), 

Améliorations Lorenz et Raz [LROl] ont proposé un algorithme polynomial pour trouver des 
bornes inférieure et supérieure dont le rapport est n. En considérant le graphe réduit d = {V, Ai) 
où Ai est l'ensemble des arcs dont le coût est parmi les i plus faibles, il suffit de trouver le i tel que 
Gi admet un chemin réalisable en consommation de ressource et Gi-i n'en admet pas. La borne 
inférieure est donc le plus grand coîit, noté C£,^+, sur Gi (tout chemin réalisable emprunte au 
moins un arc de coîit ci^f;+, sinon il existe un chemin réalisable sur Gi_i) et la borne supérieure 
est n C]^fi+ (le plus court chemin en consommation de ressource sur Gi est réalisable, de valeur 
au plus n C]^ji+). 

Cette procédure est en 0{n log^(n) + m log(n)) (en simplifiant avec Ollogm) = Oilogn)). 
De plus, les auteurs mettent en œuvre le schéma avec la valeur f{LB, UB) = (5 = 1 et p = 2. La 
complexité totale obtenue est de : 

/ m n\ 
O ym n loglog(7i) + j 

Ergun et al., [ESZ02] exploitent plus finement la qualité de ces bornes initiales pour améliorer 
encore la procédure de recherche dichotomique, en rendant dynamique la mise à jour des para- 
mètres S et B de cette procédure. Ainsi, ils démontrent que si l'on dispose en entrée de DICHO 
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Algorithme 11 : BORNE : Procédure Bornes 



Entrées : Un graphe G 

Sorties : Deux bornes LB et UB telles que LB < OPT <UB <n LB 

LR- ^ 0; 
LR+ <- m; 

tant que LR^ < LR^ — 1 faire 

LR+ + LR- 
[ 2 J' 

Calculer P^^t plus court chemin de 1 à n en consommation de ressource dans Gi = (Vf, Ai) où 
At = {(î,i) e A\c^j < Cl}; 
si (7(Pf_^,j) < b„ alors 
L ^ l; 

sinon 

L M- ^ l; 

LB +~ c^ii+ ; 
UB^C{PtS); 



U B 

de bornes LB et UB vérifiant —— ■ < n, alors, en posant : 

LB 



ô^fiLB,UB) =J^-l 



LB 



B = g{LB,UB,5) 



UB LB 

1 + 5 



La complexité de DICHO devient 0{m n) ; par conséquent, le schéma dans sa globalité se déroule 
en temps : 

Etude de la complexité Ergun et al. [ESZ02] effectuent le calcul de complexité de leur schéma 
en utilisant les valeurs exactes pour S et B. Hassin [Has92] utilise une méthode dont la complexité 
est connue pour trouver une valeur approchée des racines carrées. Ici, l'étude commence par 
choisir de bonnes valeurs approchées des paramètres du schéma suivant cette méthode avant de 
faire l'étude de la complexité. La recherche des valeurs approchées des racines carrées se fait 
comme suit : 

UB 

- trouver le premier i tel que = 2 > 

- B = LB aj_3 d'où UB^LB^ < B < UB^LB^ 
D'après le résultat de la procédure TEST, les bornes sont mises à jour : 

- TEST répond OUI : 
LB+ ^B et UB+ = UB 

^ , 1 , fUB\i UB+ fUBY 

Cela donne donc < -r- < 

\LB J - LB+ \LB) 
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TEST répond NON : 

LB+ = LB et UB+ = {l + ô) B 



Cela donne donc 



fUB\ 



UB 
< ^ < 



f 



UB 



\LB J - LB+ \LB 



Donc au pire, la plus petite diminution du rapport donne 



UB^ 

Tb^ 



< 



f UB 
\LB 



On note k le nombre 



de passages dans la boucle tant que de la procédure DICHO. Avant l'exécution de la procédure 
DICHO, le rapport du majorant sur le minorant est borné par n et après cette procédure, il est 
de 2 donc la valeur de k est bornée comme suit : 



log(7i) < log(2) 



k log 



loglog(n) < loglog(2) 



k > 



loglog(n) - loglog(2) 



log 



Donc k < 



loglog(ri) - loglog(2) 



En posant respectivement ôi, UBi et LBi les valeurs de (5, du majorant et du minorant de l'étape 
i, la complexité de la procédure DICHO s'écrit donc : 



i=l 



Oijnn) Y.^[\ 



On a aussi 



(LB, 



Donc 



< 



\UB, 



LB} 



< 



i=l 



< ir < 



UB.:" LB- 



UBfLB: 



< 



< 



UBf - LB] 
( LB, 
\UB, 



1 



LBt 



LB, 
UB, 



UB. 
1 



1 



- car pour i < A;, on a 



UB, 
LBi 



> 2 donc 



LB, 1 
1jB,,^2 



S, 



< 



i=i 



LB, 
UB, 



2 + 23 

( LB, 
\UB, 



22+23 



et on rappelle que ^ ^ < 



LB, 



UB,^i 
LB,^i 
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Finalement, cela donne 




< 2-3 car UBk < 2LB, 

3=0 

en faisant tendre k vers oo, on obtient 




La complexité de la procédure DICHO est bien en 0{m n), ce qui permet d'arriver à une com- 



5.2 Approche multicritère 

Si le cas d'une ressource se résoud efficacement par la programmation dynamique, en revanche, 
que peut-on dire du cas de deux ressources ou plus ? Bien siir, puisque décider même s'il existe 
un chemin réalisable est NP-complet, on ne peut garantir de trouver des solutions approchées 
en temps polynomial. Néanmoins, une manière de généraliser l'approche précédente consiste 
à se placer dans le cadre multicritère. il s'agit d'une relaxation du problème (on relâche les 
contraintes de consommation de ressource), mais le fait de conserver la nature multicritère (les 
consommations de ressource sont vues comme des critères à optimiser) permet peut-être de 
mieux parcourir l'ensemble des solutions ou, à défaut, d'en dessiner les contours. Dès le début 
des années 80, Hansen s'est intfessé au problème de plus court chemin bicritère en en donnant 
\m FPTAS [HanSO]. Pour des problèmes plus généraux, Papadimitriou et Yannakakis, [PYOO], 
ont proposé une démonstration géométrique de l'existence d'une frontière de Pareto approchée 
de taille polynomiale. Les auteurs ont également donné des théorèmes d'existence de procédures 
polynomialcs permettant de construire de telles frontières, qui s'appliquent notamment à notre 
problème (ces théorèmes sont d'ailleurs une forme de généralisation des résultats obtenus par 
Hassin pour le cas d'une ressource, que l'on peut considérer comme problème bicritère). Nous 
définissons dans un premier temps le problème de plus court chemin dans un cadre multicritère 
afin d'introduire la notion de frontière de Pareto approchée. Nous proposons ensuite un algorithme 
permettant de déterminer une telle frontière, puis concluons en replaçant ces résultats dans le 
cadre spécifique de la résolution de PCCC. Dans cette section, on se restreint de nouveau aux 
contraintes de ressource au puits ; de plus, les graphes sont supposés sans circuit et les données 
numériques rationnelles strictement positives. 

Plus court chemin multicritère Dans l'approche multicritère, les problèmes considérés sont 
la généralisation à plusieurs objectifs des problèmes d'optimisation classique. Pour le problème 
de plus court chemin, il s'agit donc de déterminer, dans un graphe orienté G = (V^A), un 
chemin de s h t qui optimise, non pas une fonction de coût C, mais un ensemble de fonctions 
C^, . . . , C^. Cette version du problème est notée PCC-M. Les fonctions à optimiser sont toutes 
supposées à valeur dans Q!j_. De plus, les critères peuvent être à minimiser ou à maximiser (par 
exemple, un critère souvent pris en considération pour la construction de chemins est celui de 
la longueur des chemins, que l'on cherche à maximiser). Ainsi, par la suite, on considérera R 
fonctions C^, . . . , à optimiser, sans présupposer de leur sens d'optimisation (maximiser ou 
minimiser), ni dans un premier temps de leur interprétation (critère original du problème ou 
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relaxation d'une contrainte de ressource) . Les ensembles des indices des critères à maximiser et à 
minimiser seront respectivement notes IW^'^^ et Pour tout vecteur w de M''', wi^max (resp., 

î«l^min) désigne sa restriction aux indices de Si™^^ (resp., de ^"""). 

La multiplicité des fonctions à optimiser fait qu'elles n'induisent plus a priori un ordre complet 
sur l'ensemble des solutions réalisables : (1, 2) et (2, 1) sont par exemple incomparables dans R^. 
Evidemment, il est toujours possible de définir malgré tout un ordre complet, par exemple, 
en considérant l'ordre lexicographique, ou encore, en considérant, lorsque cela est possible, une 
combinaison linéaire des critères à optimiser ; auquel cas on se ramène à un problème monocritère. 
Néanmoins, si l'on souhaite conserver la nature multicritère du problème (et c'est notre cas : la 
ressource 1 ne doit pas plus excéder sa borne que ne le doit la ressource 2), alors il faut travailler 
sur l'ordre partiel induit par la relation de dominance, déjà présentée en section 4.1. On ne 
cherche plus dès lors à déterminer une solution optimale, mais un ensemble de solutions efficaces, 
ou non dominées, appelé frontière de Pareto. 

Définition 5 Frontière de Pareto pour PCC-M ; 

Soit I = (G = (y, A); s, t e F; C\ . . . , : "P Q) une instance de PCC-M, où V désigne 
l'ensemble des chemins sur G. La frontière de Pareto de I, notée ^{I), est l'ensemble des 
solutions non dominées de I. Formellement, P{I) est l'ensemble des chemins Ps~*t qui vérifient 
pour tout chemin P' s^t S V ■ 

Sz 3rGl R ta ( e ^'"-) A (C'-(P;^,) > C^(P,^,)) 

Al ^0 , R + / (^6^"^^'^) A (C^(P,^t) > C'^(^i^t) 
Mors d« G i, . . . t.q. | ^ ^min) ^ < C'iPl^t) 



EYontière de Pareto e-approchée Potentiellement (mais pas nécessairement), que l'on se 
situe dans l'espace des solutions ou dans celui des valeurs, la frontière de Pareto est de taille 
exponentielle : déterminer cette frontière n'est donc pas envisageable. En outre, pour certains 
problèmes (et c'est notamment le cas du plus court chemin), décider même si un point donné est 
ou non Pareto optimal est déjà NP-complet. Aussi, à défaut de manipuler l'ensemble des points 
Pareto-optimaux, il pourrait être intéressant d'exhiber un sous-ensemble de taille raisonnable 
de points qui permettent de représenter, de façon approchée mais maîtrisée, tous les points de 
l'ensemble. C'est là la vocation de la frontière de Pareto e-approchée. 

Définition 6 Frontière de Pareto e-approchée pour PCC-M ; 

Soit I = (G = (y, A);s,t e F;C\ . . . : P ^ Q) une instance de PCC-M, où V désigne 
l'ensemble des chemins sur G. Un sous-ensemble S^e(X) de solutions réalisables de I est une 
frontière de Pareto e-approchée si : 

yPi^t e V, 3Ps^t e ^e{I) t.q. : 

f G'^iPi^t) < (1 + e)G''{Ps->t) sire ^""^'^ 

[ G''(P,Ui) > (1 - e)G'-{Ps^t) sire ^™ 

Cette définiton signifie que tout chemin Pg^^ a un témoin Ps~,t dans ^dl) dont le vecteur 
de peformance réalise au moins 1/(1 + e)G^ {P!._^f) si r e|^max, au plus 1/(1 — e)C^ {P!._,^) si 

r Êl^min. 

Les auteurs dans [PYOO] démontrent que tout problème multicritère admet une frontière 
de Pareto e-approchéc de taille polynomialc en |/| et en 1/e (mais exponentielle en R). Ce 
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résultat pourrait sembler surprenant, la taille d'une frontière de Pareto étant potentiellement 
exponentielle; pourtant, la preuve est relativement simple; nous la présentons ici dans le cadre 
de PCC-M. Si C^jj„ et C^^^ désignent respectivement un minorant et un majorant de C^{Ps^t) 
(on rappelle que la valeur C^{Ps^t) est supposée être rationnelle strictement positive) pour tout 
r, on découpe tout d'abord l'intervalle [Cmi„,C'^^j] en une suite ([c[,c[+i])- de intervalles, 
oià cl^i = (1 + e)c[ si r e âê"^'^^, c^^^ = (1 — e)c[ si r g pour un paramètre d'erreur e. 

Cette procédure de quadrillage de l'espace des valeurs est précisément décrite dans l'algorithme 
QUADRILLAGE. Si l'on définit la quantité MAJ comme MAJ = maxf^i{C;,„j/C;,„}, la 

complexité de cet algorithme est d'ordre O ^ ^ log{MAJ)j ^ Si la description proposée 

différiencie les paramètres d'erreur em et Sm pour les critères à maximiser et à minimiser, on 
suppose pour l'instant : ej\/ = £„i = £• 



Algorithme 12 : QUADRILLAGE : quadrillage de l'espace des valeurs 

Entrées : / instance, eA/,em £]0, 1[ erreurs pour les critères à maximiser / minimiser 

Sorties : U^^i{(cS, . • . , c^r)} discrétisation de l'espace [C^im C^naj] 

pour tous les r e faire 

pour tous les i £ {0, . . . , H^} faire 

L = (1 + SMyCrnin; 

pour tous les r £ i^™'" faire 

~ \^'~'S(l-E^) (Cmin/C'maj)! i 

pour tous les i £ {0, . . . , H^} faire 

|_ ~ (1 £m) C^j-^aj] 



Une fois le quadrillage effectué, on considère chacun des hypercubes [c^i , c^^^^] x . . . x [cfn , c^^;^] 
pour (î^, . . . , i^) e {0, . . . , iï^ — 1} X . . . X {0, . . . , H^- — 1}. Il suffit alors de choisir, lorsqu'un tel 
point existe, un point de i^e{I) P^^r hypercube. L'ensemble 3^e{I) ainsi construit est bien une 
frontière de Pareto e-approchée, de taille bornée par : 



O |^Qlog(A/AJ)^ j 



En effet, soit Pg^t un chemin de s à É ; son étiquette (C^(P,'_,j), . . . , C^(Pg^()) appartient 
à un hypercube H représenté par le point (c^i , . . . , c^) dont les coordonnées sont définies par : 

r c\r < CiP^^t) < c^^+i si r e ^"'^'^ 

1 > C^iP^^t) > c[.+i si r e ^™ 

L'hypercubc H considéré consiste alors en le produit cartésien des intervalles [c^r, c[r_|_]^[ pour 
les indices correspondant à un critère à maximiser et ]clr^i, c^r] pour les indices correspondant 
à un critère à minimiser. Il existe nécessairement une étiquette E correspondant à un chemin de 
dans H (cette étiquette pouvant être celle de P^^t) ; alors, par construction : 

< E"- et CiPl^t) < (1 + e)c\r =^ C"'(P,'_J < (1 + e)E'' si r e i^'"'^'' 

clr- > E"- et C""(P;_^J > (1 - e)clr ^ C"'(P,UJ > (1 - e)E'' si r e 
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La procédure de recherche d'une frontière de Pareto e-approchée que nous venons de discuter 
est formeUement décrite dans l'algorithme PARETO-EPSILON. Cet algorithme n'est cependant 
pas polynomial, puisqu'on ne sait a priori pas décider de l'existence d'un chemin dans un hy- 
percubc donné. Notons pour finir que cette construction reste valide si l'on choisit pour chaque 
point (c^i , . . . , cfa) un chemin P'^^i dont l'étiquette se situe dans l'hypercube plus large : 



Algorithme 13 : PARETO-EPSILON : Frontière de Pareto e-approchée 
Entrées : I instance, e G]0, 1[ erreur 
Sorties : ,^e{I) Frontière de Pareto e-approchée 
^ 0; 

U?=i{(c5,...,cïfO} ^ QUADRILLAGE(£,e); 

pour tous les (i\ . . . , i^) £ {0, . . . , iï^ - 1} x . . . x {0, . . . , //-^ - 1} faire 

si 3 un chemin Pa^t à valeur dans (^rea?"»" [c^r , c^r+i] O^g^^min [Cir_|_i, c^r] alors 



Schéma complet d'approximation de la frontière de Pareto Dans [PYOO] sont énoncés des 
théorèmes de caractérisation de la constructibilité en temps polynomial de frontières de Pareto 
£-approchées (oià, par polynomial, on entend polynomial en |/| et 1/e, exponentiel en R) ; ces 
constructions sont qualifiées de schémas complets d'approximation pour le problème multicritèrc. 
La version multicritèrc d'un problème d'optimisation admet un schéma complet d'approximation 
de sa frontière de Pareto ssi l'on peut décider en temps polynomial la version décision approchée 
(de type GAP) du problème multicritèrc. En particulier, si les fonctions à optimiser sont linéaires, 
discrètes et que la version monocritère exacte est décidable en temps pseudopolynomial, alors 
le problème multicritèrc admet un schéma complet d'approximation. S'agissant du problème de 
plus court chemin, on sait par la programmation dynamique décider en temps pseudopolynomial, 
pour une instance / et une valeur B, s'il existe sur / un chemin de valeur B ; on en déduit donc que 
PCC-M admet un schéma complet d'approximation de sa frontière de Pareto. Nous proposons ici 
un tel schéma, qui est une forme de généralisation au cadre multicritèrc des algorithmes proposés 
dans la section précédente. 

Principe. Dans le cadre monocritère à une ressource, pour une instance / et une borne B données, 
on ne sait pas décider en temps polynomial s'il existe un chemin de coiit au plus B tout en 
consommant au plus 6„ quantité de ressource. En revanche, on peut pour une erreur s donnée 
décider en temps polynomial (en / et en 1/e) s'il existe un chemin de coût au plus (1 + e)B et de 
consommation au plus 6„, ou si tout chemin, soit coûte au moins B, soit consomme strictement 
plus que b„ . Cette notion de test approché se généralise naturellement à plusieurs fonctions : on 
ne sait pas, pour un point (c^i , . . . , c^) donné, décider s'il existe un chemin Ps^t dont l'étiquette 
{C^{Ps-.t), ■ • ■ , C^iPs^t)) vérifie pour tout critère CiPs^t) > si r e CiPs^t) < 

si r € i^™™, ou si tout chemin n'atteint pas ou excède strictement l'une au moins des bornes c[r. 
En revanche, on peut décider en temps polynomial pour des erreurs Em et £,„ données s'il existe un 
chemin respectant ces bornes, ou si tout chemin est en-dcça de {l + SM)clr pour un certain critère 
à maximiser ou au-delà de (1 — Sm)clr pour un certain critère à minimiser. De nouveau, ce test 
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approché est obtenu en travaillant sur une instance / modifiée : on ramène les données numériques 
de l'instance / initiale à un ordre polynomial par une procédure d'échelonnage et d'arrondis 
dépendant du point (c^i , . . . , c[r) et des degrés de précision £m, £m ; au prix de l'erreur introduite, 
on gagne le déroulement polynomial sur / des algorithmes exacts qui sont de complexité théorique 
pseudopolynomiale. 

Algorithme de décision exacte. Nous présentons ici l'algorithme ACYCLIQUE-M qui est une 
adaptation au cadre multicritère de l'algorithme à fixation d'étiquette 5. Cet algorithme renvoie, 
pour une instance / de PCC-M et un vecteur B de Q^, l'ensemble (a des étiquettes E G non 
dominées vérifiant : 

ACYCLIQUE-M renvoie également les chemins associés aux étiquettes de S', en ne conservant 
toutefois qu'un chemin par étiquette. Selon l'analyse faite en section 4.2, la complexité de cet 
algorithme est d'ordre : 

n I i-1 \ 

iî + ^ ^ NjR + C%M + NnR 

i=2 \j=l j 

Or, si les données sont entières, si M désigne un majorant du nombre de valeurs possibles 
(considérer par exemple, si l'on note Cmaj = max^^ilC'maj} 6t Cmin — niiii^=i{C'min}j ^ = 
(n — l){Crnaj — Cmin + 1) < 0{nCmaj))^ le nombre Nj d'étiquettes en tout sommet j est borné 
par M^. Si l'on ne conserve qu'un chemin par étiquette en chaque sommet, la complexité de 
ACYCLIQUE-M est alors d'ordre au plus : 

O {RM"" {j? + nAf«)) = O {Rn^''+' Ci^,^) 

L'algorithme ACYCLIQUE-M permet ainsi de décider en temps pseudopolynomial, pour une 
instance / de PPC-M aux données numériques entières et un vecteur B, s'il existe un chemin 
dans l'hyperespace (g)rg^max[i?'', -|-oo[®re5?inin[0, B^], ou si tout chemin excède B^ pour un certain 
critère r S âê"^"^, ou n'atteint pas B'' pour un certain critère r G âê"^'^^. 



Procédures d'échelonnage et d'arrondis et de test approchée. Pour se ramener d'un ordre pseu- 
dopolynomial à un ordre polynomial, on échelonne les données numériques de l'instance / par 
l'algorithme SCALING-M dont la complexité est d'ordre : 

O (mR log (^)) = O (n^R log (^)) 

L'arrondi effectué pour obtenir des données entières sur / a pour conséquence que la décision n'est 
plus exacte, mais seulement approchée, sur l'instance initiale (algorithme TEST-M). Considérons 
l'instance / obtenue à partir de / après l'appel à la procédure SCALING-M pour les paramètres 
B = (Cj\ , . . . , c^), £m et Ejn- On voit facilement que l'on a les relations suivantes entre les 
étiquettes sur les instances I et I d'un chemin Ps_»t : 

C^iP.^t) > Bm Vr G J"»^'^ r C{Ps^t)\3ê^.. > B|^-x 
C'-{Ps^t) < Bm Vr e 1 C(P,_t)|5^„.„ < B|5^.„ax 

3r e C"^(Ps^t) < Bm j 3r e C'-(P,^t) < (1 + eM)B- 

3r e ^-i", C^iPs^t) >Bm ^ 1 3r e C^{P,^t) > (1 - em)B^ 
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Algorithme 14 : ACYCLIQUE-M : Algorithme pour les graphes acycliques 

Entrées : I = {G = {V, A); c\ . . . , : 4 ^ Q) instance de PPC-M, B = (B\ . . . , B^) vecteur borne 
Sorties : Ensemble des étiquettes {E^ , . . . ,E^) non dominées vérifiant < ~ ^ ^ ^min 

// ETIQ(i) est l'ensemble des étiquettes du sommet i 

ETIQil) ^ (0) 

pour tous les i allant de 2 à n — 1 faire 

ETIQii) ^ 0; 

pour tous les j £ predeœsseur{i) faire 
pour tous les i? e ETIQ{j) faire 
si Vr G + c'ji < alors 

E' ^{E^ + c^^,,Wre{l,...,R}); 
ETIQii) ^ ETIQii) U {£"}; 

ETIQii) ^ ParetoiETIQii)); 
ETIQin) ^ 0; 

pour tous les j G predecesseurin) faire 
pour tous les Ê G ETIQij) faire 

si Vr G + c;„ < et Vr G + 4„ > B" alors 

£' ^(iî'- + c,-„,Vr G {!,..., J?}); 
ETIQin) ^ ETIQin) U {£'}; 

ETIQin) ^ ParetoiETIQin)); 
Retourner ETIQin) ; 



Algorithme 15 : SCALING-M : Procédure d'échelonnage et d'arrondis multicritère 



Entrées : / = (G = (V, A); , . . . , : A - 
Sorties : liB ,eM , Em) instance de PPC-M 
pour tous les les (i, j) G A faire 

r— 1- 



instance de PPC-M, B G (Q+)-^, SM.Em G]0, 1[ erreurs 



Ba 



pour tous les r G . 



EA/B'- 



faire 

Bm 



pour tous les r G . 



£mB 



faire 
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Par ailleurs, puisque les valeurs numériques sur / sont bornées par max{_BM, Bm} = 0{n/ min{ejv/, 
le déroulement de la procédure ACYCLIQUE-M sur / pour le paramètre B défini par i?|^max = 
Bm et iJi^min = Bjn devient polynomial en n et l/e^/, l/e™, d'ordre : 



O 



min{£A/,e„}2^ 



En conséquence de ces deux dernières observations, on peut décider, pour une instance / 
d'ordre n, un point {c\i , . . . , c^) et deux paramètres d'erreur £M,£m, en temps polynomial en n 
et en max{l/eA/, l/e™} : 

(Si 3Ps^f d'étiquette dans l'hyperespace (K),-g^n,ax [c[r, +oo[(8irGâ?n""[0, c[r] 
1 Ou si l'hyperespace [(1 + f M)c-'r , +oo[(g)rg<2?mi„[0, (1 - est vide. 



Algorithme 16 : TEST-M : Procédure de test approchée sur PCC-M 
Entrées : / instance, B = {B^ , . . . , B^) vecteur, £M,£m G]0, 1[ erreurs 
Sorties : 

• OUI et un chemin Ps^t venhant < ^,r ,n \ ^ r,r vj ^ ^min 

[ C (-Ts-^t) < B Vr € £ig 

• NON si tout chemin Ps-»t vérifie C^{Pa-,t) < (1 + eM)B^ pour un certain critère r à maximiser ou 
C'^{Ps~,t) > (1 — em)B^ pour un certain critère r à minimiser. 

Î{B, eM,em) ^SCALING-M(/, B, em, e^); 
Bm 



n 




n 




; Bm ^ 




£m 





B' ^ {B"- = Bm, Vr e B"- = B^, & .^'"'"); 

ETIQ ^ACYCLIQUE-M(/, B'); 
si ETIQ / alors 
|_ Retourner OUI et un chemin d'étiquette dans ETIQ 

sinon 

L Retourner NON 



Schéma complet d'approximation Le schéma conçu consiste à opérer un quadrillage de l'espace 
des valeurs, puis à appeler la procédure de test approché en le coin de chaque hypercube issu de ce 
quadrillage. Précisément, l'algorithme PARETO-PCC-M solhcite la procédure QUADRILLAGE 
pour les erreurs Em = vT+e — 1 et = 1 — yî— £ (on ne fait plus le quadrillage sur £, puisque 
le test en chaque hypercube n'est pas exact, mais approché). Il appelle ensuite la procédure 
TEST-M sur chaque coin (c^i , . . . , c^), en invoquant les mêmes paramètres d'erreur em et e„i- 
Les chemins trouvés par les appels successifs à TEST-M constituent la frontière de Pareto £- 
approchée. 

Soit Pg^t un chemin de s à t et soit (C^(Pg'_^j), . . . , {P'g_^^)) son étiquette, on considère 
le point B = (c^^i , . . . , c^) caractérisé par : 

< C'-(Pi_t) < (1 + EaiK- si r e ^'"^'^ 
> CiP^^t) > (1 - em)<- si r e 
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Algorithme 17 : PARETO-PCC-M : Frontière de Pareto e-approchée en temps polynomial 



Entrées : / instance, e £]0, 1[ erreur 

Sorties : !3^e{I) Frontière de Pareto e-approciiée 

eu ^ y/l + £ - 1 ; 
Sm ^ 1 - VI - e; 

U?=i{(c5,...,cïf.)} QUADRILLAGE(/, em, e™); 

pour tous les (i\ . . . , i^) £ {0, . . . , Jï^ - 1} x . . . x {0, . . . , //^ - 1} faire 

B^{c]^,...,c%)- 

si TEST-M(I ,B,eM ,£m) renvoie un chemin Ps^t alors 
L ^.(/) ^ ^.(/)U{P.^t}; 



Si TEST-M(/, i?, £jv/, £m) a renvoyé un chemin Ps-»t, alors ce chemin vérifie : 

r ^^(Ps^*) > ^ C"'(P,U) < (1 + EmX,. < (1 + e)C^iPs^t) si r e .^^^^ 
\ C'\Ps^t) < => C-{P's^t) > (1 - EmX^ > (1 - e)C'\Ps^t) si r e ^-"^ 

Sinon, l'existence de P'g_,i assure que TEST-M a nécessairement renvoyé un chemin Pg^t 
pour le point i?_ défini par : 

B!:=c[._iVr = l,...,iî 
Ce chemin vérifie bien de nouveau : 
r C-{Ps^t) > T+^<. ^-(n^J < (1 + SM^C-iPs^t) = (1 + e)C'-(n^t) si r e 
\ C'^(P.^t) < j^cl. => C^{P's-.t) > (1 - e™)'C'-(P,^t) = (1 - e)C''{Pg^t) si r e 

La complexité de l'algorithme PARETO-PCC-M est dominée par celle des appels successifs 
à TEST-M. Le nombre de ces appels est de l'ordre de O ( ((1/e) log(MAJ))^ ) (considérer que 



\/l + £ — 1 ^ 1/26, 1 — \/T~~ë ~ l/2£). Or, la complexité de TEST-M somme la complexité 
temporelle de SCALING-M appfiquée à / et celle de ACYCLIQUE-M appliquée à /. On obtient 
au total une complexité d'ordre : 

' '\og{MAJ)\" f ^_fn^ , Rn^^+^\\ _ ^ f f\og{MAJ)\'^ Rn^^+^\ 



£ / \ V £ / £^^ J I \\ £ J I 

Frontière de Pareto et PCCC La vision multicritcre du problème de plus court chemin 
contraint à R ressources peut s'énoncer comme le problème consistant à trouver un chemin 
Ps-,t qui minimise chaque composante du vecteur (^C{Ps^t),T^{Ps^t), ■ ■ ■ ,T^{Ps^t)) ■ Si l'on 
dispose d'une frontière de Pareto £-approchée cela signifie que tout chemin P^^f de s à 
t est représenté par un chemin Ps^t de la frontière qui consomme au plus 1/(1 — £) la 
consommation de P^^t sur chaque critère. En particulier, si Pg_,f est réalisable, on déduit que 
Ps_»t est 1/(1 — £)-réahsable. 

6 Synthèse 

Le tableau 1 résume les complexités théoriques des algorithmes de plus court chemin les plus 
connus ainsi que des algorithmes exacts ou approchés pour le PCCC vus dans ce rapport. 
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Tab. 1. Différents algorittimes et leur complexité 



Nombre 
de 


lïypothèse 
sur 

ic gl cLJJlic 


Auteur 
ou 

TVT/^wn ^ 1 ' Q 1 rv/^T"i ^" Vi tvi û 
INUlli U.t; i digUi ILlllllt: 


Complexité en temps 





C^J G Q 


Bellman[Bel58] 


0{mn) 





G acycliquc 


Bellman 


0{m + n) 





Cij £ Q+ 


Dijkstra[Dij59] 


0{{m + n) log(n)) 





c,j e Q+ 


Dijkstra (tas de Fibonacci) 


0{m + n log(n)) 


f 


Cij e 


Programmation dynamique [DB03] 


0{mbi) 


ReNt 


c^j e Q 


Programmation dynamique 


0(e") 




C^J G Q 


k plus courts chemins [Epp98] 
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Le but de cette étude est d'étudier la possibilité de générer des chemins contraints en temps 
maîtrisé, dans l'optique d'accélérer certains schémas de génération de colonnes impliquant le 
PCCC comme sous- problème. Deux avantages pourraient être tirés d'une telle démarche : d'une 
part, être en mesure de générer une population diversifiée de chemins lors de la phase d'initialisa- 
tion; d'autre part, accélérer les premières itérations du schéma. Les méthodes les plus répandues 
à ce jour pour résoudre le PCCC, notamment dans ce cadre spécifique, sont celles de type pro- 
grammation dynamique, qui offrent des solutions exactes, en temps et en espace exponentiel (ce 
sont encore ces mêmes méthodes qui ont inspiré les algorithmes approchés) . Si nous avons évoqué 
quelques méthodes de simplification des instances (section 3), celles-ci ne garantissent pas pour 
autant la diminution de la complexité théorique. Quelle alternative avons-nous donc, quitte à 
s'affranchir de l'optimalité. voire, de la réalisabilitc ? 

La recherche de plus courts chemins dans le cadre spécifique du déroulement d'un schéma de 
génération de colonnes pose deux problèmes de nature différente. D'une part, celle de la difficulté 
intrinsèque du problème contraint, à partir de 2 ressources : si l'on souhaite rester dans un ordre 
polynomial de complexité, il faut s'affranchir de l'exactitude en termes d'optimalité et en termes 
de réalisabilité. D'autre part, les coûts manipulés sur les arcs ne sont plus les coiits initiaux de 
l'instance, mais des coûts réduits issus de la résolution du (PMR), qui peuvent être négatifs : les 
hypothèses faites sur la structure de coût de l'instance initiale ne tiennent donc plus au cours 
des itérations. Or, ces hypothèses interviennent conjointement dans les calculs de complexité et 
des facteurs d'approximation. De sorte à exploiter les résultats d'approximation, il faudrait donc 
être en mesure de les étendre à un cadre plus général (notamment, coûts rationnels négatifs et 
positifs) ; en outre, il faudrait utiliser, voire, définir, un cadre pertinent pour l'approximation, 
l'approximation classique supposant des valeurs positives. 

La diversification des chemins (notamment, dans la phase d'initialisation) sous-entend un 
certain balayage de l'espace des solutions. Parmi les approches que nous avons présentées, la 
recherche des k plus courts chemins d'une part (section 4.3), d'une frontière de Pareto e-approchée 
d'autre part (section 5.2), vont dans le sens d'une telle exploration. Dans le premier cas, on cherche 
intensivement les meilleurs chemins contraints en terme de coût, tandis que dans le second cas, 
on balaye sporadiquement tout l'espace. 

Pour un nombre de ressources R donné, l'algorithme d'Eppstein permet de trier en temps 
0{m, + Rn log(ri,)) tous les chemins du graphe suivant le coût (initial, ou réduits dans le cadre 
de la génération de colonnes). Une fois le tas constitué, il suffit pour trouver le chemin contraint 
optimal de dépiler les chemins jusqu'à trouver un chemin réalisable (ou, dans le cadre d'un 
schéma de génération de colonnes, quand le critère de coût choisi pour la construction du tas 
est le coût réduit, l'optimalité est prouvée s'il n'existe pas de chemin réalisable de coût négatif). 
Rien n'assure néanmoins que les k plus courts chemins, pour k = 0{p{n)) où p polynôme, 
contiennent ne serait-ce qu'un chemin réalisable pour les ressources (ni un chemin de coût positif 
pour démontrer l'optimalité du schéma). Par ailleurs, en triant sur le seul critère de coût, on risque 
de générer des chemins peu diversifiés. Ainsi, en se limitant à chercher un nombre polynomial de 
chemins, on ne garantit pas l'obtention de chemins pertinents ; en prenant au contraire k = 0{e"), 
la construction du tas reste polynomiale en temps et en espace, et l'énumération des chemins 
devient exponentielle en temps mais pas espace : on n'obtient alors une alternative à la résolution 
exacte par programmation dynamique qui est exponentielle en temps et en espace. 

A partir de 2 ressources, l'obtention de solutions réalisables n'étant plus garantie, l'approche 
multicritère semble prometteuse. Bien plus que la recherche de la frontière de Pareto e-approchée 
elle-même, c'est le quadrillage de l'espace des valeurs et l'utilisation de la programmation dyna- 
mique sur les instances modifiées qu'il faudra chercher à exploiter : balayer l'espace des valeurs 
possibles et chercher une solution dans chaque hypcrcube semble être une bonne stratégie en vue 
de la diversification des chemins dès l'initialisation. 
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